home *** CD-ROM | disk | FTP | other *** search
/ The X-Philes (2nd Revision) / The X-Philes Number 1 (1995).iso / xphiles / hp48hor2 / leastsq.doc < prev    next >
Text File  |  1995-03-31  |  3KB  |  75 lines
  1. LEASTSQ - Least squares curve fitting 
  2. From: dnnant01@ucthpx.uct.ac.za (Anthony K. Donno) 
  3. Newsgroups: comp.sources.hp48 
  4. Date: 23 Sep 91 03:27:56 GMT 
  5.  
  6. The feature that I use most in the HP48 stats menu is that of curve fitting, 
  7. however if you have data that does not fit the standard curve types offered 
  8. by the calculator, a good fit cannot be obtained. 
  9.  
  10. The following program fits a polynomial equation to a set of data points in a 
  11. similar fashion to the stats menu. It takes a matrix off the stack that has at 
  12. least two columns in it (this is not checked), allows you to specify which 
  13. column is x data and which is y data, and then fits a curve to the data using 
  14. a polynomial of an order that you specify. Please note that more than 2 
  15. columns are allowed, but the program only works with the two that are 
  16. selected. 
  17.  
  18. The program uses a standard least squares fit, but operates fairly quickly 
  19. (4.5 mins to fit a 15th order to 21 data pairs). Note that an option allows 
  20. you to specify whether the program must plot the fitted curve, or just return 
  21. the equation. The program does not recalculate the polynomial, unless you 
  22. change the x or y columns or the order of the polynomial required. 
  23.  
  24. [If you want an exact polyfit (rather than least-squares approximation), see 
  25.  POLYFIT on this disk.  -jkh-] 
  26.  
  27. The following menu is used: 
  28.  
  29. XCOL   - takes an integer off the stack and makes that column the x data 
  30.          column 
  31. YCOL   - takes an integer off the stack and makes that column the y data 
  32.          column 
  33. ORDER  - takes an integer off the stack and makes it the order of the poly- 
  34.          nomial that is required 
  35. PLOT   - calculates (if necessary) and plots the fitted polynomial 
  36. FIT    - calculates (if necessary) and returns the polynomial that fits the 
  37.          data 
  38. EXIT   - returns back to the LEAST directory CST menu (you can put whatever 
  39.          you want in this menu although a suggested menu is given) 
  40.  
  41. The following is a list of the programs that make up LEAST: 
  42.  
  43. LEAST  - main program (menus ...) 
  44. PLOT   - plotting function 
  45. FIT    - curve fitter (brains) 
  46. POLY   - converts a list of numbers on the stack into a algebraic polynomial 
  47.          of the order given in level one (the numbers are the coefficients) 
  48. SETRNG - sets the range for the plot 
  49. XYR    - redraws the top line of the LEAST screen (XCOL: 1 YCOL: 2 ORDER: 3) 
  50. DELEQ  - deletes the current equation EQ 
  51.  
  52. The following variables are produced by LEAST: 
  53.  
  54. ORDR  - order of the polynomial 
  55. YC    - y data column 
  56. XC    - x data column 
  57. dDAT  - summation data (same as the stats menu) 
  58. EQ    - current equation 
  59.  
  60. Checksum : # 63743d 
  61. Bytes    : 1950.5 
  62.  
  63. I hope that this program is of use to others, and would appreciatiate it if 
  64. you would mail any comments to me. 
  65.  
  66. ============================================================================== 
  67. Anthony Donno                          | 
  68. Department of Electrical Engineering   |   Caesar entered on his head 
  69. University of Cape Town                |   he had a helmet on each foot 
  70. New South Africa                       |   he had a sandal in his hand 
  71.                                        |   he had his trusty sword to boot 
  72. e-mail : dnnant01@ucthpx.uct.ac.za     | 
  73. phone  : +27 +21 6892420               | 
  74. ============================================================================== 
  75.